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1.  Introduction 


C«d«a 

- FTi  aad/or 

Special 


Diet 


The  suboptimal  Hoo  problem  can  been  solved  via  the  the  solution  of  two  algebraic 


Riccati  equations  (Doyle  et  al.,  1989).  In  the  finite  horizon  case,  expressions  for  a  subopti¬ 
mal  controller  have  also  been  derived  in  the  time-varying  setting  (Subrahmanyam  1992a), 
without  resorting  to  various  transformations,  such  as  the  ones  given  in  (Safonov  et  al., 
1989).  In  the  time-invariant  case,  the  solutions  of  the  dynamic  Riccati  equations  eventu¬ 
ally  converge  to  constant  matrices.  Also,  efficient  algorithms  for  the  computation  of  the 
infimal  ifoo  norm  have  been  given  in  the  finite  horizon  case  (Subrahmanyam  1992b, c,d). 
Our  computational  experience  indicates  that  for  time-invariant  systems,  the  infimal  ffoo 
norm  can  be  very  nearly  approached  using  full  state  feedback,  whereas  with  output  feed¬ 
back,  the  suboptimal  value  that  gives  a  viable  controller  is  generally  much  higher  than  the 
infimal  value  (Subrahmanyam  1992e). 

It  is  difficult  to  design  a  suboptimal  controller  taking  into  account  parameter  varia¬ 
tions  in  the  system  matrices.  It  is  doubtful  whether  a  noniterative  solution  exists  to  the 
performeince  robustness  problem.  The  reason  for  this  is  the  fact  that  when  one  designs 
a  suboptimal  Hoo  controller  at  the  nominal  values  of  the  parameters,  the  performance  of 
the  controller  under  a  variation  in  the  system  matrices  is  a  priori  unknown.  In  fact,  in 
our  treatment  here,  knowledge  of  the  controller  is  essential  for  the  computation  of  the 
variation  of  performance  imder  variations  in  the  system  matrices.  It  is  tacitly  assumed 
that  the  controller  and  observer  matrices  designed  at  a  nominal  point  are  fixed  even  under 
parameter  uncertainties  since  we  have  no  way  to  measure  the  variations  in  the  system 
matrices  in  general. 
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To  solve  the  performance  robustness  problem,  we  need  to  have  an  idea  of  the  variations 
in  performance.  To  this  end,  we  derive  an  expression  for  the  variation  in  performance  for 
a  specific  controller  in  terms  of  variations  in  the  system  matrices.  If  the  performance 
degradation  for  the  range  of  parameter  variations  is  unsatisfactory,  we  need  to  redesign 
the  suboptimal  controller. 

The  technical  note  is  organized  as  follows.  In  Section  2  we  formulate  the  finite  horizon 
Hoo  performeuice  robustness  problem.  Section  3  presents  a  siunmary  of  the  equations 
involved  in  the  design  of  a  suboptimal  H^o  controller.  Section  4  presents  a  procedure 
for  the  computation  of  actual  performance  of  the  suboptimal  controller.  In  Section  5, 
a  formula  for  the  degradation  in  performance  owing  to  parameter  variations  is  given. 
Section  6  presents  an  iterative  procedure  for  the  design  of  a  suboptimal  controller  which 
has  adequate  performance  robustness.  An  example  of  the  control  system  design  which 
illustrates  the  usefulness  of  the  theory  is  given  in  Section  7. 


2.  Problem  Formulation 

Let  the  n-dimensional  time- varying  system  be  given  by 

X  =  A{t)x  +  Bi(t)u  +  B2{t)v^  x(<o)  =  0, 

z  =  C(<)z  -I-  D(t)u  +  E{t)v, 
y  =  C2(f)i  +  D2{t)u  -f 


(2.1) 

(2.2) 

(2.3) 

(2.4) 
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where  R  and  W  are  assumed  to  be  positive  definite  and  the  superscript  *  denotes  a  matrix 
or  vector  transpose.  Computational  techniques  for  the  evaluation  of  Xopt  are  given  in  the 
papers  by  Subrahmanyam  ( 1992b, c,d).  The  suboptimal  Hoo  problem  can  be  stated  as 
follows.  Given  A  <  Xopti  find  full  state  and  output  feedback  controllers,  if  these  exist,  for 
which 


•  /o  \v*Rvdt 

rmn-^^V? - > 

/o 


A. 


(2.5) 


The  above  problem  has  been  solved  in  Subrahmanyam  (1992a)  and  a  summary  of  the 
design  equations  will  be  given  in  Section  3. 

The  suboptimal  performance  robustness  problem  can  be  stated  as  follows.  Given 
A  <  Xopt,  find  state  and  output  feedback  controllers,  if  these  exist,  for  which 


J/'  \v‘Ib>dt 

/(T 


mm 


>A, 


even  under  specified  variations  in  the  system  matrices  in  (2.1)-(2.3). 


3,  Feedback  Solutions 


In  this  section,  we  give  equations  for  a  suboptimal  Hoo  design. 

derivation,  see  Subrahmanaym  (1992a,1991). 

The  system  is  given  by 

For  the  details  of 

X  =  A(i)x  +  Bi{t)u  +  B2{i)Vf  x(0)  =  0, 

(3.1) 

z  =  C{t)x  +  D(t)u  +  E{t)vy 

(3.2) 

y  =  C2{t)x  +  2?2(f)u  +  E2{t)v. 

(3.3) 
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Let 


X  •  So 

Xopi  =  maxmm  -%r-7 - , 

'  «  •’^0  J^^z*Wzdt 

and  let  A  <  Xopt-  Also,  let  •  •  ■ , be  defined  by 


(3.4) 


z'Wz  =  x'Wix  +  Tx’Wiu  +  u’WiU  +  2i"H'4»  +  v'WfV  + 


(3.5) 


Let 

Zi  =  +  E^ii  +  E*w,  (3-6) 

and  let  ,  •  •  • ,  be  defined  by 

z\R~^z\  =  e*W^ic  +  2e*W^2U  +  +  2c*H^4t5  +  ty'TVstij  +  2ti*U^6«>-  (3.7) 

We  make  the  following  assvimptions. 

(a)  Let  the  final  time  T  satisfy  the  assumption  in  Lemma  4.2  of  Subrahmanyam  1992a. 

(b)  Assume  that  for  each  t  €  [0,r],  R  —  AWj  is  positive  definite  and  W3  +  AW6(il  - 
XW^y^W;  is  invertible. 

(c)  Also,  assume  that  for  each  t  e  (0,r],  W~^  —  XW^  is  positive  definite  and  W3  + 
XWeiW-^  -  AW5)"^W6  is  invertible. 

The  relevant  controller  equations  are 


Q  =  {R-XW3)-\ 

(3.8) 

A  =  {W3  +  AW«ftWV}-', 

(3.9) 

Ui  =  A{B;  +  AWeftBJ), 

(3.10) 

t/2  =  -A(w;  +  AWenw;), 

(3.11) 

V^=XQi-B;  +  W:Ui), 

(3.12) 

V2  =  Afi(W7  +  w;u2), 

(3.13) 
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P  +  P{A  +  B1U2  +  B2V2)  +  (A*  -  W2VX  -  WMP 

+P{BxUi  +  B2Vx)P  -  (Wx  +  W2U2  +  W4V2)  =  0,  P(T)  =  0.  (3.14) 

Note  that  the  above  equation  is  symmetric. 

The  relevant  observer  equations  are 

$  =  -  AW5)-^  (3.15) 

r  =  {W^3  +  XWe^W^}-\  (3.16) 

i  =  A  +  +  (3-17) 

B  =  ~(UxP  +  U2),  (3.18) 

C  =  C2  +  E2ViP  +  E2V2,  (3.19) 

51  =  r(C  +  XW^^DB),  (3.20) 

52  =  -T{W;  +  XWe^W:),  (3.21) 

Ti  =  X^{-DB  +  W;Sx),  (3.22) 

T2  =  X^W:  +  W^S2).  (3.23) 

Y  =  (A-  W2SX  -  W4Ti)Y  +  Y{A*  +  C'Sa  +  B’PTj) 

+y(C*Si  +  B*VTi)Y  -  {Wx  +  1^252  +  W4T2),  y(0)  =  0.  (3.24) 

Note  that  the  above  equation  is  symmetric. 

A  suboptimal  controller  is  given  by 

q  =  Aq  +  Bx{U\Pq  +  U2q)  +  B2(VxPq  +  V29)  +  B{Cq  +  D2^  ~  y),  (3.25) 

Z  =  (Siy  +  52)*,  (3-26) 

u  =  {UiP  +  U2)q.  (3-27) 
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In  the  case  of  time-invariant  systems,  the  solutions  of  the  Riccati  equations  (3.14) 
and  (3.24)  eventually  converge  to  constant  matrices.  Note  that  the  full  state  feedback 
controller  is  simply  given  by  u  =  {UiP  +  U^)!. 


4.  Computation  of  Performance 

We  consider  the  output  feedback  case.  The  full  state  feedback  case  is  covered  by  the 
ensuing  analysis  as  well.  The  closed  loop  system  is  given  by 


(i\_(  A  B,{UiP  +  U2) 

\qj  ~  \ -LCi  U  -  LDiiUiP 


?«)(:)*(-&)•■  <•» 


where 


U  =  A  +  BiiUiP-\-U2)  +  BiiViP  +  Vi)  +  Lie  +  D^iUiP  +  U2)). 


(4.2) 


Note  that  all  the  matrices  can  be  time-varying.  We  now  specialize  to  the  case  where 


x(0)  =  g(0)  =  0. 


(4.3) 


Let 


For  the  controller 


the  performance  is  given  by 


z  —  Cx  +  DiUiP  -J- 1^2)9  "b  Ev. 


u  =  iU^P  +  U2)q, 


(4.4) 


(4.5) 


\v*Rvdt 

min  - .  (4.6) 

^z*Wzdt 

The  above  value  is  strictly  greater  than  A  by  our  design  procedure.  For  time-invariant  sys¬ 
tems,  it  is  convenient  to  consider  constant  solutions  of  the  two  dynamic  Riccati  equations 
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and  take  the  matrices  U\P  and  S\Y  +  5^  to  be  fixed.  In  some  cases,  it  is  possible 
that  the  performance  may  be  below  A  because  of  the  disregard  for  the  time  dependence. 
Nevertheless,  in  such  cases,  the  performance  can  be  artificially  improved  by  designing  the 
controUer  at  a  higher  value  of  A  than  the  one  spedfied. 

We  now  describe  a  procedure  to  compute  the  value  of  (4.6).  Let  (4.1)  be  rewritten  as 

X,  =  A,x»  +  x,(0)  =  0,  (4.7) 


where  v  needs  to  be  chosen  to  minimize 

_ \v*Rvdt _ 

/o^{  +  x*Q2V  +  |v*Q3u}  dt 

Note  that  the  denominator  of  (4.6)  can  put  in  the  form  of  the  denominator  of  (4.8)  by 
virtue  of  (4.4). 

This  problem  has  previoxisly  been  solved  by  Subrahmanyam  (1992f,  1990)  and  con¬ 
ditions  guaranteeing  the  existence  of  a  minimizing  v  are  also  reported  in  these  references. 
For  the  sake  of  completeness,  we  present  here  the  detedls  concerning  the  characterization  of 
a  minimizing  v  and  the  computation  of  the  performance  corresponding  to  the  minimizing 

V. 

the  following  theorem,  we  give  conditions  that  are  satisfied  by  an  optimal  v  which 
minimizes  (4.8)  subject  to  (4.7). 


Theorem  4.1.  Consider  the  system  given  by  (4.7)  and  (4.8).  Let 


A  =  min 


J^^v*Rvdt 


fo  +  ^0Q2V  +  jVQsv)  dt 
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Assume  that  (R(t)  -  ^Qzit))  *  exists  for  all  t  &  [0,T].  If  (x,,t;)  is  optimal,  then  there 
exists  a  nonzero  p(t)  such  that 


^  =  -a:p-Xqiz.-xq,v,  p(r)  =  o, 


and 


v(t)  =  (R-  XQsr'fxCnx.  +  B.V)- 

Proof.  If  V  minimizes  (4.8),  then  it  also  minimizes 

rT  ,  ,T 


(4.10) 


(4.11) 


/  Iv’Rvdt-xj  {~x:Q,x,  +  z:QiV+\v"Q3v}dt.  (4.12) 

The  theorem  now  follows  from  the  maximum  principle  (Pontryagin  et  al.,  1962).  □ 

Let 


and 


A  =  A,  +  XB.{R-XQ,)-'Ql 

b  =  b.(r-xq,)-'b:, 

C  =  -XQj  -  X^Q,(R  -  XQ3)-‘Q;. 


(4.13) 

(4.14) 

(4.15) 


The  variables  satisfy  a  two  point  boimdary  value  problem  given  by 

(7) -a 


with 


1.(0)  =  0,  p(T)  =  0. 


(4.16) 


(4.17) 


A  A 

We  now  give  a  criterion  for  the  estimation  of  A.  Notice  that  A  gives  a  measrure  of 
performance  of  the  ifoo  optimeJ  controller  under  worst-case  conditions. 
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Theorem  4.2.  Let  A  be  the  smallest  positive  value  for  which  the  boundary  value  problem 
given  by  (4.16)  and  (4.17)  has  a  solution  (x„p)  with  jJ^{^XtQi^»d-XgQ2V+^v*Q3v}  dt  > 
0,  where  v  =  (R  —  AQs)”^  +  13*p}.  Then  A  is  the  minimum  value  of  (4.8),  Xg  is  an 

optimal  trajectory  and  v  =  (jR  —  AQa)”^  +  B*p}  is  an  open  loop  optimal  exogenous 

input. 


Proof.  It  is  clear  from  Theorem  4.1  that  if  t;(<)  minimizes  (4.8),  then  it  satisfies  ^4.16) 
and  (4.17),  with  A  being  the  minimum  value  of  (4.8).  Now  suppose  (x„p)  ^  0  satisfies 
(4.16)  and  (4.17)  for  some  A.  Let  v  =  {R—\Q3)~^  {XQ2Xa+B*p].  In  the  following  equations 
(  ,  )  denotes  an  inner  product. 

We  have 


r  {{R-XQ3)v,v)dt=  rC\Q*2X„v)dt+  r{B:p,v)dt.  (4.18) 

Jo  Jo  Jo 


By  equation  (4.7),  the  second  integral  of  (4.18)  can  be  written  as 

f  {B*p,v)dt=  f  {p,B,v)dt=  f  (p,x,  -  Agx,)dt.  (4.19) 

Jo  Jo  Jo 

An  integration  by  parts  and  equations  (4.10)  and  (4.17)  yield 

/  {p,Xg  -  AgXg)dt  =  X  {Qix,,x,)dt  + X  {x,,Q2v)dt.  (4.20) 

Jo  Jo  Jo 

Substituting  (4.20)  in  (4.18),  we  get 

I  v*Rvdt  =  xf  {x*QiX, +  2x*Q2V  +  v*Q3v}  dt.  (4.21) 

Jo  Jo 

Thus,  the  cost  associated  with  t;  is  A  provided  that  the  right  side  of  (4.21)  is  nonzero. 

Hence,  if  (x,,p)  is  a  nontrivial  solution  of  the  boundary  value  problem  given  by  (4.16)  and 
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(4.17)  for  the  smallest  parameter  A  >  0  with  Jq{xIQiX,  +  2xlQ2V  +  dt  >  0,  then 

A  is  the  optimal  value  and  (x«,p)  is  an  optimal  pair.  □ 

Note  that  the  boundary  value  problem  (4.16)-(4.17)  has  a  solution  with  a  nonvanishing 
denominator  for  (4.8)  for  at  most  a  coimtably  infinite  values  of  A.  Theorem  4.2  gives  a 
sufficient  condition  for  an  open  loop  exogenous  input  to  be  optimal.  Theorems  4.1  and  4.2 
completely  characterize  the  open  loop  worst  case  exogenous  input. 

Making  use  of  the  transition  matrix,  the  solution  of  (4.16)  may  be  expressed  as 


( ^-(0^  _  f  ^n(<,0)  $i2(f,0)\  f  x,(0)\ 

V p(t) ) - v^2i(<.o)  $22(<,o); \m)' 


(4.22) 


Equation  (4.17)  yields 


$i2(T,0)/>(0)  =  x,(r),  (4.23) 

^22(r,0)p(0)  =  0.  (4.24) 

In  view  of  (4.24)  and  (4.16)-(4.17),  we  have  det($22(r, 0))  =  0  if  and  only  if  the  solution 
(x,,p)  of  (4.16)-(4.17)  is  not  identically  zero.  Thus,  we  need  the  least  positive  A  which 
makes  det($22(T',  0))  =  0  and  the  denominator  of  (4.8)  positive.  This  can  be  very  easily 
obtained  by  doing  a  search  with  A  over  an  interval  on  which  there  is  a  change  in  the  sign 
of  the  determinant. 

We  foimd  the  following  algorithm  to  be  ntunerically  more  stable  since  numbers  of 
lesser  magnitude  are  involved  in  the  computation  of  the  transition  matrices  in  (4.25).  We 
have 

(  kt/  r/2)«(r/2, 0)  (  )  .  (4.25) 
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«-'(T,T/2)  =  I*), 

and 

Making  use  of  x,(0)  =  p{T)  =  0,  we  have 


I/12P(0)  = 
I/22P(0)  = 

Thus 

det 

V  1/22 


inx.iT), 

(4.26) 

6ix,(T). 

(4.27) 

o 

II 

(4.28) 

Thus  we  need  the  least  positive  A  which  malces  the  above  determinant  zero. 


5.  Performance  Variation 

In  this  section,  we  develop  a  formula  for  the  variation  of  A  when  there  are  parameter 
variations  in  the  system  matrices.  The  system  is  given  by 


with  V  chosen  to  minimize 


it  =  AtXt  +  BfV,  x,(0)  =  0, 


Jo  \v*Rvdi 


+  \v*Qzv] 

From  the  theory  of  Section  4,  we  have  a  boundary  value  problem  given  by 

i,  =  Az,  +  Bp, 
p  w  Car,  -  A*p, 


(6.1) 


(5.2) 


(5.3) 

(5.4) 
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with 

x.(0)  =  0,  p{T)  =  0.  (5.5) 

Note  that  A  is  the  minimal  value  of  (5.2).  Because  of  variations  in  the  matrices 

of  the  original  system  (2.1)-(2.3),  there  will  be  corresponding  variations  in  the  matrices 

A(t),  jB(t),  and  C{i).  Let  the  elemental  dependent  variations  in  A,B,  and  C  be  denoted 
by  &4,  SB,  and  SC  respectively.  For  simplicity  of  notation,  we  will  derive  the  variation  in 
performance  in  terms  of  6A,6B,  and  SC.  We  denote  by  Si  the  variation  in  A  owing  to  the 
parsimeter  variations. 

Let 

r  =  {R-XQ,)-'.  (5.6) 

We  have 

SA  =  1i+SJi,  (5.7) 

SB  =  /2  +  SXJ2,  (^-S) 

SC  =  h+Sj3,  (5.9) 

where 


7i  =  M.  +  xsb,tq;  +  xB,rs  r  -  a«?3)tq;  +  xb^tsq;,  (5.10) 

Ji  =  -XB.TQsTQ;  +  B.TQl  (5.11) 

I2  =  sBsTb;  +  B,r  SR-  a«?3)tb;  +  B,rsB:,  (5.12) 

J2  =  -B.TQsTB;,  (5.13) 

/3  =  -XSQi  -  A2«?2Tg;  -  X^Q2rS  R  -  XSQ2)rQl  -  >?Q2’t8Ql,  (5.14) 

h  =  -Qx  -  2XQ2'tQl  +  A"g2TQ3Tg;.  (5.15) 
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Let  Xi  and  pi  represent  the  variation  in  x,  and  p  due  to  SA,6B,  and  SC.  FVom  (5.3)- 
(5.5),  we  have  the  following  equations  that  are  satisfied  by  Xi  and  pi. 


ii  =  Axi  -f  Bpi  -f-  (/i  -l-  £XJi)i  -|-  {I2  +  £XJ2)p,  (5.16) 

Pi  =  Cxi  —  A* Pi  +  (I3  -f  AJa)!  —  (Ii  d-  SXji)*p,  (5.17) 

a:i(0)  =  0,  p,(r)  =  0.  (5.18) 

Theorem  5.1.  The  variation  S\  in  performance  is  given  by 

S\  =  ~  Jo  ^  dt  ^ 

Jq  {^x*Qix,  +  x*Q2V  +  \v*Q3v}  dt' 

where  v  =  T{5*/)  +  XQ^x). 


Proof.  From  (5.17)  we  obtain 

rT  rT  rT 


f  x*pi  dt  =  f  x*Cxi  dt-  f  x*A*pi  dt 
Jo  Jo  Jo 

+  f  x:il3  +  ^J3)x,dt-  f  x:{Ii+^Jiypdt.  (5.20) 
Jo  Jo 

Integrating  the  left  side  of  (5.20)  by  parts  and  making  use  of  (5.3)  and  (5.5),  we  obtain 

—  /  p*iBpdt=  I  x*Cxidt+  I  x^Iz  +  ^J3)x,dt  —  I  x*{Ii  +  ^JiYpdt.  (5.21) 
Jo  Jo  Jo  Jo 

By  (5.4),  the  first  integral  on  the  right  side  of  (5.21)  is  written  as 

/  x:Cxidt=  /  {p  +  A*pyxidt.  (5.22) 

Jo  Jo 

An  integration  by  parts  and  equations  (5.16)  and  (5.18)  3rield 
fT  _  fT  fT  ^  fT 

I  x*Cxidt  =  —  I  p*Bpidt—  I  p*(Ii  +  ^Ji)xgdi  —  I  p*{l2  +  ^J2)x:tdt.  (5.23) 
Jo  Jo  Jo  Jo 
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Substituting  (5.23)  in  (5.21)  and  simplifying,  we  obtain 

&  =  (5  24) 

/o  {2x:  JiV  +  p*J2p  -  x;  J31,}  dt 

It  cem  be  easily  shown  (see  Lemma  1  of  Subrahmanyam  1992f,  for  example)  that  the 

fry 

denominator  of  (5.24)  is  equal  to  Jjj  {x*QiXt  +  2xlQ2V  +  dt.  □ 


6.  Performance  Robustness  Problem  Solution 

We  now  give  an  iterative  procedure  to  solve  the  performance  robustness  problem  posed 


in  Section  2.  The  equations  are  given  by 

X  =  A{t)x  +  Bi{t)u  +  B2(t)v,  x(0)  =  0, 

(6.1) 

z  =  C(t)x  +  D(f)u  +  E{t)v, 

(6.2) 

y  =  C2(t)x  +  D2{t)u  +  E2{t)v, 

(6.3) 

and 

X  .  So  \v*Rvdt 

Aop(  =  max  mm  ^  - . 

«  L;^*Wzdt 

(6.4) 

The  suboptimal  performance  robustness  problem  is  to  find  a  controller, 

if  one  exists, 

for  which 

•  So  k^*Rvdt  -  - 

mm  -^r-f - >  A,  A  <  A„pf, 

/o  i^*Wzdt 

(6.5) 

under  variations  in  the  system  matrices  involved  in  (6.1)-(6.3). 
The  iterative  steps  are  as  follows. 


1.  Set  A  =  Ai,  A,'  <  Aopt,  and  design  a  suboptimal  Hgo  controller  based  on  the  design 
equations  of  Section  3. 


14 


NAWCADWAR-TN-92041-60 


2.  For  this  controller,  let 


A  =  min 


fJ'^v*Rvdt 


fo  +  ^iQ2V  +  iv*Q3v}  dt 


(6.6) 


This  can  be  computed  using  the  results  of  Section  4.  Note  that  A  >  A.  If  A  <  A, 
increase  the  value  of  A  and  go  through  Steps  1  and  2  until  A  >  A. 

3.  For  the  range  of  allowable  variations  of  the  matrices  involved  in  (6.1)-(6.3),  find  the 
worst  value  of  6A  using  (5.19).  Note  that  (5.19)  is  linear  in  the  parameter  variations. 
Thus,  let  6  =  min  fix.  The  controller  of  Step  1  solves  the  performance  robustness 
problem  in  case  A  +  fi  >  A.  If  this  is  the  case,  stop  the  iteration  and  use  the  controller 
of  Step  1. 

4.  If  A  4-  fi  <  A,  choose  a  suitable  Aj+i,  Aj+i  <  Xopt,  and  set  A  =  A,+i.  Go  back  to  Step  1. 
Note  that  A, 4.1  may  be  bigger  or  smaller  than  A,-  depending  on  the  specific  problem 
being  solved. 

Alternately  we  can  get  A  and  fi  =  minfiA  for  a  range  of  values  of  A  and  then  pick 
a  value  of  A  that  solves  the  performance  robustness  problem.  In  case  we  cannot  find  a 
value  of  A  which  satisfies  the  performance  robustness  requirement,  then  the  requirement 
is  too  stringent  and  the  performance  robustness  problem  does  not  have  a  solution.  The 
requirement  can  be  relaxed  either  by  lowering  A  or  by  suitably  shrinking  the  range  of 
allowable  parameter  variations. 


7.  An  Example 

We  now  illustrate  the  theory  with  an  example.  The  system  equations  are  given  by 
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(6.1)-(6.4)  with 

^=(J  _\)'  =  o)- 

^={1 9’  -^=(0  0)- 

C2  =  (i  0),  i>2  =  o,  f;2  =  (o  1),  R  =  100/2,  w  =  i2, 

where  I2  is  the  2-dimensionaJ  identity  matrix.  Let 


Qi  =&4(1,1), 


02  =  1), 


03  =  682(2, 1). 


Assume  that  a,  €  [— 0.1,0.1],i  =  1,2,3.  Given  a  A,  the  problem  is  to  find  a  controller  for 
which 


r„^v*Rvdi 

mm-:^ - > 

o/o  f^z*lVzdt 


(7.1) 


under  the  above  variations  in  A(l,  l),Ri(l,  1)  and  82(2, 1). 

We  can  find  Ao^*  given  by  (2.4)  using  the  technique  given  in  Subrahmanyam  1992b,c,d. 
For  this  example,  ^opt  =  261.4016.  A  suboptimal  full  state  feedback  controller  can  be 
designed  choosing  A  arbitrarily  close  to  Aopt-  However,  to  obtain  a  viable  suboptimal  output 
feedback  design,  A  needs  to  be  chosen  much  lower  than  Xopt-  Let  us  form  Table  1  which 
gives  the  performance  for  various  values  of  A  and  the  variations. 
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TABLE  ] 

:  Results  of  the  Performance  Robustness  Example 

A 

A 

Ml 

M4 

A-^ 

1 

5.4601 

-0.7034 

-3.0474 

-0.0428 

2.3868 

2.4127 

2 

5.8783 

-0.6757 

-3.2158 

0.0169 

2.5232 

2.6625 

3 

6.3303 

-0.6428 

-3.3935 

0.0837 

2.6670 

2.9368 

4 

6.8179 

-0.6045 

-3.5808 

0.1577 

2.8186 

3.2371 

5 

7.3431 

-0.5605 

-3.7778 

0.2392 

2.9781 

3.5653 

6 

7.9079 

-0.5107 

-3.9847 

0.3285 

3.1455 

3.9232 

7 

8.5140 

-0.4549 

-4.2010 

0.4254 

3.3207 

4.3130 

8 

9.1630 

-0.3933 

-4.4254 

0.5296 

3.5026 

4.7376 

9 

9.8562 

-0.3258 

-4.6543 

0.6398 

3.6886 

5.2019 

10 

10.5944 

-0.2534 

-4.8779 

0.7528 

3.8717 

5.7165 

We  now  explain  the  procedure  to  obtain  a  particular  row  of  the  table.  For  example, 
for  A  =  3,  we  find  the  output  feedback  controller  using  the  equations  given  in  Section 
3.  Note  that  we  ignore  the  time  dependence  of  the  solutions  (3.14)  and  (3.24)  and  only 
consider  the  constant  solutions.  For  A  =  3, 


UiP  +  Ui  =  i  -2.4249  - 1 .0062 ) , 


ViP  + Vj  = 


L  = 


The  closed  loop  system  (4.1)  is  given  by 


given  by 

-2.4249 

-1.0062 

0 

0 

-4.1773 

-0.0062 

-0.2745 

-0.9849 

(7.2) 


v.’ith  x(0)  —  9(0)  =  0.  Going  through  the  procedure  given  in  Section  4,  the  actual  perfor¬ 
mance  A  of  the  output  feedback  controller 


u  =  {UiP  +  U2)q 


(7.3) 
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A  A 

c£ui  be  obtained  as  the  first  positive  value  of  A  which  satisfies  (4.28).  This  value  is  A  = 
6.3303.  Using  (4.26)-(4.27)  and  selecting,  for  example,  the  first  component  of  />(0)  as  1,  we 
can  obtmn  p(0).  Thus,  we  can  solve  the  initial  value  problem  (5.3)-(5.4)  with  r»(0)  =  0 
and  the  value  of  p(0)  as  obtained  above. 

Assuming  that  the  accuracy  of  linear  expressions  is  acceptable,  03,03)  which  is 

the  variation  in  performance  can  be  obtained  by  using  (5.19).  Since  (5.19)  is  linear  in  the 
p2irameter  variations,  the  worst  degradation  in  the  value  of  A  is  at  one  of  the  8  vertices 
of  (01,03,03).  Because  of  the  linearity  of  AX(oi, 03,03)  and  the  symmetric  nature  of  the 
variations,  we  need  only  compute  at  four  vertices.  In  Table  1, 

6Xi  =a(0.1,0.1,0.1), 

6X3  =£X(0.1,-0.1,0.1), 

^3  =M(0.1,0.1,-0.1), 

6X4  =  6A(-0.1,0.1,0.1). 

Also  6  =  min,  ±^Ai.  Thus  A  —  6  in  Table  1  gives  the  greatest  lower  boimd  of  the  actual 
performance  of  the  controller  under  parameter  variations.  This  value  is  2.9368. 

Going  back  to  (7.1),  if  A  =  3,  looking  at  Table  1,  A  can  be  4  or  higher.  If  A  =  4,  the 
controller  design  can  be  performed  with  A  =  7,  or  with  a  higher  value  of  A.  If  A  =  5,  we 
can  select  A  at  9  or  higher. 


8.  Conclusions 

In  this  technical  note,  an  approach  is  presented  for  the  solution  of  the  finite  horizon 
Hoo  performeince  robustness  problem  imder  parameter  variations.  A  linear  expression  for 
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the  degradation  of  performance  is  given  in  terms  of  variations  in  the  system  matrices.  An 
exBimple  which  illustrates  the  methodology  is  also  given. 
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